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ABSTRACT 

As part of the Greenland Sea Project Woods Hole Oceanographic 
Institution and Scripps Institute of Oceanography deployed six 
acoustic tomography transceiver moorings to measure variability of 
the Greenland Sea gyre through a cooling cycle from September 1988 
to August 1989. Using a set of Greenland Sea acoustic tomography 
data provided by Woods Hole Oceanographic Institution this thesis 
investigated the importance of incorporating acoustic path changes 
in the construction of the tomographic inverse solution. 

A comparison of the inverse solutions for changes in sound 
speed using non-corrected and corrected acoustic multipaths was 
conducted. Although the two inverse solutions are qualitatively 
Similar, significant quantitative differences exist. These 
differences indicate that it is necessary to account for changes in 
the acoustic multipaths for the generation of accurate Greenland 
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I. INTRODUCTION 


A. OCEAN ACOUSTIC TOMOGRAPHY 

The purpose of tomography is to remotely sense the 
interior structure of a desired medium. This may be done 
through a variety of methods. X-rayS are used in medical 
computer assisted tomography (CAT) scans and seismic waves are 
introduced by surface sources in the study of the interior 
structure of the earth. Ocean acoustic tomography is 
accomplished by transmitting sound waves through the sea, a 
medium which is relatively transparent to sound. 

Munk and Wunsch (1979) first introduced the procedure of 
using acoustic inverse techniques to monitor ocean basins for 
mesoscale fluctuations. The application of acoustic 
tomography to monitor the interior structure of the ocean has 
two parts. The first part of acoustic tomography is to solve 
the "forward" problem. Sound propagates through the ocean 
from a source to a receiver along the acoustic multipaths. 
Measured travel time changes of the sound signals along these 
paths constitute the data of this underwater acoustic 
observing system. Given a reference sound speed field, it is 
necessary to predict the geometry and travel time of sound 
waves along these acoustic multipaths. This can be done by 


using ray tracing computer programs. The second part of 


acoustic tomography is solving the “inverse" problem. Travel 
time changes of sound signals are used to infer changes in the 
sound speed field. The inversion of travel time changes can 
be accomplished by a linear, minimum mean square error 
estimator as has been shown by Cornuelle et al. (1985), Chiu 


et al. (1987), and Joseph (1991). 


B. GREENLAND SEA PROJECT 

The Greenland Sea Project (GSP) is a five year (1987- 
1992), international scientific study to understand the large 
scale, long term interactions of air, sea, and ice in the 
Greenland Sea. The five key elements of GSP are to study the 
air-sea-ice interaction, ocean ventilation, ocean circulation 
and mixing, atmospheric energetics, and biological processes. 

As part of GSP, an array of six tomography transceivers 
were deployed in the Greenland Sea in the fall of 1988, by 
Woods Hole Oceanographic Institution (WHOI) and Scripps 
Institute of Oceanography (SIO). The role of ocean acoustic 
tomography was to assist in the studies of ocean ventilation 
and circulation. The array of transceiver moorings is 
designed to measure changes in the integrated properties of 
the Greenland Sea gyre through a winter cooling season 
(Greenland Sea Science Planning Group, 1986). Figure 1 shows 
the deployment site and geometry of the array. Each mooring 
transmitted an acoustic signal six times a day every third 


day, during the period of September 1988-August 1989. Table 


1 provides the transceiver locations. The latitude and 
longitude of the moorings were converted to an XY plane by 
Joseph (1991). This operation was conducted by placing a 240 
km by 240 km grid over the six transceivers. The grid was 
centered on Mooring 6. Figure 2 depicts the grid system used 


to convert the mooring position to XY coordinates. 


C. OBJECTIVE 

The objective of this thesis is to investigate the 
importance of incorporating acoustic path variability in the 
tomographic inverse solution. Travel time of sound signal 
along an acoustic path is affected by changes in the sound 
speed along the path. In addition, it is also affected by 
changes in the path. In the construction of the inverse 
solution, the path is usually held stationary and only the 
changes in sound speed are considered. It is the intent of 
this study to determine if path geometry changes need to be 
considered in the inversion of the Greenland Sea tomography 


data. 


D. APPROACH 

The thesis objective was met by comparison of the inverse 
solutions for changes in sound speed obtained using non- 
corrected and corrected multipaths. Only the transmissions 
between Moorings 4 and 5 were used in my investigation. The 


acoustic travel time data were provided by WHOI. The 
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Figure 1. Deployment Site of the Greenland Sea Tomography Array. (from 
Greenland Sea Science Planning Group, 1986) 
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Table I: TRANSCEIVER LOCATIONS. (from Joseph, 1991) 
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"forward" problem was solved by using the Multiple-Profile Ray 
Tracing Program (MPP), a ray theory based algorithm, to 
calculate the acoustic multipaths. MPP was developed by 
Spofford (1973). The output from MPP was then interfaced with 
the inversion code to calculate the estimated change in sound 
speed (6c). The inversion code was developed by Chiu (1987) 
and modified by Joseph (1991). 

All the inversions performed in this thesis used a twelve 
day average of the observed data. That is, twelve-day- 
averaged travel time changes were used as input to the 
inversion code. The first method of data inversion used non- 
corrected multipaths. The non-corrected multipaths were 
calculated by MPP using a range independent Sound Velocity 


Profile (SVP). This SVP was measured by WHOI and SIO during 


oa 
& 
&» 


oe 
~ 
oe 


—_ 
E 
a4 
oo” 
v& 
oD 
c 
cS 
fees 


06 120 144 
range (km) 





Figure 2. PLAN VIEW OF ARRAY GEOMETRY. (from Joseph, 1991) 


the deployment of the moorings. The multipaths were held 
constant through a period from day 270 to day 342. The study 
will focus on this period. This procedure will be referred to 
as the static ray approximation. The second method of data 
inversion used corrected (i.e., updated) multipaths. The 
initial multipaths were exactly the same as those used in the 
static ray approximation. However, as time progressed, the 
background profiles were updated. This created a set of range 
dependent SVPs for MPP to calculate updated multipaths. This 
procedure will be referred to as the updated ray geometry. If 


the 6c maps generated by the two methods show significant 


differences, then it would be necessary to account for changes 
in the acoustic multipaths in the inversions of the Greenland 


Sea tomography data. 


II. DESCRIPTION OF APPROACH 

This chapter details the systematic approach used to 
achieve the thesis objective as stated in the introduction. 
MPP was used to predict the acoustic multipaths or eigenrays. 
All eigenrays were calculated from a source at Mooring 4 to a 
receiver at Mooring 5. Only this mooring pair was used in 
this study. Eigenrays are the paths that acoustic signals 
travel for a particular source-receiver geometry. In the 
Greenland Sea, eigenrays are usually grouped in bundles of 
four. The initial angles of the rays ina bundle are within 
a few tenths of a degree from each other. Two of the rays in 
a bundle start in a downward direction and the other two start 
in a upward direction. Upon arrival at the receiver two rays 
in the bundle arrive from beneath the receiver and two from 
above. For our inverse analysis we have used the four rays in 
the 13 degree bundle, containing three lower turning points, 
the 10 degree bundle, containing four lower turning points, 


and the slowest arrival from the zero degree path. 


A. INTERFACE MPP TO INVERSION CODE 

MPP was used to calculate eigenrays within an initial 
angles of +20°. Eigenrays were traced from the source to the 
receiver in both range independent and range dependent modes. 


Multiple SVP profiles at a designated interval of 24 km 


between profiles were used in each range dependent run. MPP 
linearly interpolated at 30 evenly spaced intervals between 
the designated profiles. For all runs a flat bottom at a 
depth of 5000 m was assumed. No bottom-interacting eigenrays 
were used for the subsequent inverse analysis. 

MPP’s output provides the initial angle (degrees), the 
number of discretized ray coordinates (integer variable), and 
travel time (s) of each eigenray. The output also provides at 
unevenly spaced horizontal range intervals, the ray 
coordinates (range (km) versus depth (km)), grazing angle 
(radians), and sound speed (km/s) along each eigenray. 

The output data from MPP was rearranged into an acceptable 
format for the inversion code. The rearrangement used cubic 
splines to interpolate the ray coordinates from a variable 
number of unevenly spaced horizontal range intervals into 500 
evenly spaced intervals. 

The inversion code was set up by Joseph (1991) to map 
sound speed perturbation 6c(x,y,z) 1n a 240 km x 240 km x 3 km 
ocean volume centered on Mooring 6. The horizontal dimensions 
were divided into ten 24 km segments. The vertical dimension 
was divided into five layers of variable thickness. The 
layers were centered at 100 m, 200 m, 450 m, 750 m, and 1500 
m depths, respectively. These divisions and layering 
discretized the ocean volume into 500 three dimensional boxes. 


It is essential for the estimator to accurately determine the 


amount of path length a ray has in each box based on the input 
ray coordinates. 

To incorporate path dynamics in the inversions of the 
observed travel time changes (12 day averages), the estimated 
change in sound speed (6c) calculated at a 12 day time step 
was used to update the background profiles and hence the ray 
geometry for the inversion at the next time step. This 
procedure will be discussed in greater detail later in this 


chapter. 


B. ESTIMATING ABSOLUTE POSITION ERROR 

It 1S important to estimate the absolute range error so 
that the corresponding travel time bais con be removed from 
the data. To determine the absolute range error between the 
source and the receiver the differences in the mean travel 
time of the observed signals and the travel time predicted by 
MPP for the 13 degree bundle was used. The observed mean 
travel times (T,) were obtained by averaging the data over a 
period of 11 months. The predicted travel times (T.) were 
Calculated using the background profile shown in Figure 3, in 
a range independent mode. 

To correct the horizontal range between the source and 


receiver, the following relation was used 


8x=o(T,-T)) (1) 
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where 6x is the horizontal range error and c is a mean sound 
speed. Multiple MPP runs were made varying &x until the 
calculated and observed travel times matched each other for 


the 13 degree bundle. 


C. TRAVEL TIME TIME SERIES 

The data file provided by WHOI was used to construct a 
travel time time series for this study. Each transmission day 
consisted of six transmissions. These gave six travel times 
for each individual eigenray. The travel times were corrected 
for relative mooring motion and averaged over each 
transmission day. They were then averaged over twelve day 
intervals. The last day of each interval was used as the 
first day of the next twelve day period. These twelve day 
travel time averages were used for both the static ray 
approximation and updated ray geometry inversions. 

Throughout the remainder of this thesis, day XXX, such as 
day 270, will be used to denote a twelve day period centered 


at that day. 


D. STATIC RAY APPROXIMATION 

There are two primary ways in which the travel time of 
sound signals are affected. The first is by fluctuations in 
the sound speed and the second is by fluctuations in the ray 
paths. The neglect of the second effect in the inversions of 


Gata over the period of the experiment constitutes the static 
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Figure 3. BACKGROUND SVP AND CALCULATED EIGENRAYS USED IN THE 
STATIC RAY APPROXIMATION INVERSION. 


ray approximation. 

Nine eigenrays were calculated by MPP using the range 
independent background profile shown in Figure 3. The 
Calculated eigenrays are also presented in this figure. The 
rays were held constant for each of the inversions over the 


period considered. 


E. UPDATED RAY GEOMETRY 
To account for the fluctuations in ray paths over the 
experimental period, I updated the eigenrays after the 
inversion of the data at each time step in the time series. 
The initial inversion at day 270 utilized the same range 


independent profile and eigenray geometry as displayed in 


12 


Figure 3. An estimated éc map was generated by the inversion 
code. This map was used to update the range independent 
background profile. A series of range dependent SVPs resulted 
from these corrections. Updated ray paths were then 
calculated by MPP using the range dependent profiles, and were 
used for the next inversion. In updating the background and 
ray geometry, 6c was linearly interpolated between boxes and 
eigenrays were retraced. This process was repeated for each 


inversion in the time series. 
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III. INVERSION RESULTS AND DISCUSSION 


A. ABSOLUTE POSITION ERROR 

The 13 degree ray bundle was used for the range 
correction. This bundle spent the least amount of time in the 
near surface layer which contained the greatest variability 
and thus would be least susceptible to its effects. Table III 
shows that the travel times calculated by MPP based on the 
range given by Joseph (1991) were on the average 219 ms faster 
than the observed travel times. An initial estimate of 320 m 
for the error in horizontal range was calculated. Several MPP 
runs were conducted in which the position of the source was 
fixed and the horizontal range to the receiver was increased 
until the optimal range correction was obtained. The final 
results listed in Table III demonstrate that with a correction 
of +322 m MPP was able to calculate travel times of the 13 
degree bundle to within an average of 1.98 ms of the observed 
mean travel time. The addition of this +322 m to the 
horizontal range would ensure that the differences between the 
observed and calculated travel time was due mainly to ocean 


variability. 


B. INVERSION USING STATIC RAY APPROXIMATION 
Inversions using the static ray approximation were 


conducted using the procedures previously outlined. Figure 3 
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Table H: TRAVEL TIME BIAS AND CORRECTION 
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contains the range independent background profile used by MPP 
to calculate the associated eigenrays. Figure 4 is a map of 
the estimated éc from the inversion for day 270. This map is 
the same for both inversion procedures considered in this 
study. This figure demonstrated the greatest amount of 
Variability takes place in the upper 1000 m. The sound speed 
minimum occurs at approximately 100 m. The source and receiver 


are located within the vicinity of the sound speed minimun. 
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Figure 4. ESTIMATED éc FOR DAY 270 USING THE STATIC RAY 
APPROXIMATION. (contour interval is 1 m/s) 
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Figure 5. ESTIMATED éc FOR DAY 282 USING THE STATIC RAY 
APPROXIMATION. (contour interval in 1 m/s) 
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Figure 6. ESTIMATED éc FOR DAY 294 USING THE STATIC RAY 
APPROXIMATION. (contour interval is 1 m/s) 
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Figure 7. ESTIMATED éc FOR DAY 306 USING THE STATIC RAY 
APPROXIMATION. (contour interval is 1 m/s) 
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Figure 8. ESTIMATED éc FOR DAY 318 USING THE STATIC RAY 
APPROXIMATION. (contour interval is 1 m/s) 
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Figure 9. ESTIMATED éc FOR DAY 330 USING THE STATIC RAY 
APPROXIMATION. (contour interval is 1 m/s) 
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Figure 10. ESTIMATED éc FOR DAY 342 USING THE STATIC RAY 
APPROXIMATION. (contour interval is 1 m/s) 


The éc maps estimated using the static ray approximation, 
shown in Figures 4-10, indicate an overall cooling of the 
ocean volume in the time period from day 270 through day 342. 
The upper 1000 m contains the largest amount of variability 
over this period. The ocean below 1500 m remains relatively 
stable. 

There is a positive Sc cell centered at approximately 84 
km in the horizontal. This cell increases in intensity 
between days 270 and 294. The cell decreases in intensity 
after day 294. This decrease is continuous through day 342. 
A second positive éc cell is centered at approximately 140 km 
in the horizontal. This cell decreases in intensity between 


270 and 294. Between days 294 and 342 the intensity remains 


9 


relatively constant and the horizontal area covered spreads 
out. There is a negative 6c cell centered at approximately 
168 km in the horizontal. The intensity of this cell appears 
to fluctuate throughout the entire study period. Finally, two 
negative éc cells are located at approximately 80 and 144 km 
in the horizontal and between 500 and 2000 m depths. The 
intensity of these cells remain constant from day 270 through 


day 342. 


C. INVERSION USING UPDATED RAY GEOMETRY 


Inversions using updated ray geometry were conducted using 
the procedures previously outlined. Figure 3 shows the 
initial range independent background profile and _ the 
associated eigenrays used. The odd numbered figures between 
Figures 11-22 are the updated range dependent SVPs and the 
associated eigenrays calculated by MPP. The even numbered 
figures between Figures 11-22 are the estimated 6c maps 
generated directly by the inversion code. These éc estimates 
were used to update the range dependent background profiles at 
the respective time steps. In order to facilitate the 
comparison with the static ray inversion results we show in 
Figures 23-28 the corresponding cumulative éc maps. The 6éc 
displayed in these latter maps are the changes with respect to 


the initial range independent background profile. 
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Figure 11. BACKGROUND SVP AND UPDATED EIGENRAYS FOR DAY 282. 
Starting in the upper left and moving counterclockwise; SVPs at 0 km, 24 km, 48 km, 
72 km, 96 km, and 123 km. 
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Figure 12. ESTIMATED 6c USED TO UPDATE EIGENRAYS FOR DAY 282. 
(contour interval is 1 m/s) 
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Figure 13. BACKGROUND SVP AND UPDATED EIGENRAYS FOR DAY 294. 
Starting in the upper left and moving counterclockwise; SVPs at 0 km, 24 km, 48 km, 
72 km, 96 km, and 123 km. 
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Figure 14. ESTIMATED 6c USED TO UPDATE EIGENRAYS FOR DAY 294. 
(contour interval is 1 m/s) 
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Figure 15. BACKGROUND SVP AND UPDATED EIGENRAYS FOR DAY 306. 
Starting in the upper left and moving counterclockwise; SVPs at 0 km, 24 km, 48 km, 
72 km, 96 km, and 123 km. 
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Figure 16. ESTIMATED éc USED TO UPDATE EIGENRAYS FOR DAY 306. 
(contour interval is 1 m/s) 
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Figure 17. BACKGROUND SVP AND UPDATED EIGENRAYS FOR DAY 318. 
Starting in the upper left and moving counterclockwise; SVPs at 0 km, 24 km, 48 km, 
72 km, 96 km, and 123 km. 
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Figure 18. ESTIMATED 6c USED TO UPDATE EIGENRAYS FOR DAY 318. 
(contour interval is 1 m/s) 


24 


‘ 
-1 
-1 
-1 
-f 
-? 
°2 
-2 
°? 
-2. 


ee ¢e 8 
a w»~ @eeexLewegeeeeiw”n @ 


GErTyu ite) 
i Ie Ae a ae Pe ae 
eeeenwwesek 





Figure 19. BACKGROUND § SVP AND UPDATED EIGENRAYS FOR DAY 330. 


Starting in the upper left and moving counterclockwise; SVPs at 0 km, 24 km, 48 km, 
72 km, 96 km, and 123 km. 
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Figure 20. ESTIMATED 6c USED TO UPDATE EIGENRAYS FOR DAY 330. 
(contour interval is 1 m/s) 
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Figure 21. BACKGROUND SVP AND UPDATED EIGENRAYS FOR DAY 342. 


Starting in the upper left and moving counterclockwise; SVPs at 0 km, 24 km, 48 km, 
72 km, 96 km, and 123 km. 
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Figure 22. ESTIMATED éc USED TO UPDATE EIGENRAYS FOR DAY 342. 
(contour interval 1s 1 m/s) 
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Figure 23. ESTIMATED CUMULATIVE éc FOR DAY 282 USING THE 
UPDATED RAY GEOMETRY. (contour interval is 1 m/s) 
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Figure 24. ESTIMATED CUMULATIVE éc FOR DAY 294 USING THE 
UPDATED RAY GEOMETRY. (contour interval is 1 m/s) 
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Figure 25. ESTIMATED CUMULATIVE 6c FOR DAY 306 USING THE 
UPDATED RAY GEOMETRY. (contour interval is 1 m/s) 
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Figure 26. ESTIMATED CUMULATIVE éc FOR DAY 318 USING THE 
UPDATED RAY GEOMETRY. (contour interval is 1 m/s) 
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Figure 27. ESTIMATED CUMULATIVE 6c FOR DAY 330 USING THE 
UPDATED RAY GEOMETRY. (contour interval is 1 m/s) 
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Figure 28. ESTIMATED CUMULATIVE éc FOR DAY 342 USING THE 
UPDATED RAY GEOMETRY. (contour interval is 1 m/s) 
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The 6c maps, Figures 23-28, estimated using the updated 
ray geometry indicate that the largest amount of variability 
is now confined within the upper 2000 m of the ocean volume. 
The region below 2000 m has only slight changes. 

Comparison of the 6c maps generated using the updated ray 
geometry with those generated using the static ray 
approximation indicate that the two results are qualitatively 
similar. However, there are quantitative differences. The 
positive and negative éc cells correspond to approximately the 
same locations in the two inversions. Between days 270 and 
306 the features appear to be approximately the same size in 
both the horizontal and the vertical directions. However, the 
features appear to be stronger in the inversion using the 
updated ray geometry. Between days 306 and 342 the features 
appear to be both larger and stronger in the inversion using 
the updated ray geometry. 

An analysis of the eigenrays indicate that the 13 degree 
ray bundle remains stable over the period considered in this 
study. These rays spend the least amount of time in the area 
of high variability. The 10 degree ray bundle appears to be 
more variable and the zero degree ray is the most variable 
one. The variability of the eigenrays appears to increase as 
the amount of time the ray spends in the upper portion of the 


ocean volume increases. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

The MPP raytraces indicate the presence of an absolute 
range error between the source at Mooring 4 and the receiver 
at Mooring 5 of approximately +322 m. Accounting for the 
absolute position error 1s important in reducing the amount of 
error introduced into the inverse solution. 

If fluctuations in the ray path were having negligible 
effects, Sc estimates generated by the two methods (i.e., 
static ray approximation and updated ray geometry) should be 
approximately the same. Although the two inverse solutions 
are qualitatively Similar, Significant quantitative 
differences exist in the Sc estimate generated by the two 
inversion methods. 

It was also noted that the eigenrays were sensitive to the 
strong near-surface sound speed variability. The effects on 
the elgenrays were more pronounced the longer a ray spent in 
the near surface layer. 

My conclusion is that it is necessary to account for the 
variability of the acoustic multipaths in the inversion of the 
Greenland Sea tomography data. This can be accomplished by 
updating the background profiles and ray geometry after every 
transmission day. Doing this will provide more accurate maps 
of the Greenland Sea gyre circulation through a winter cooling 


season. 
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